Specific heat of a one-dimensional interacting Fermi system: the role of anomalies 
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We re-visit the issue of the temperature dependence of the specific heat C(T) for interacting 
fermions in ID. The charge component CdT) scales hnearly with T, but the spin component CsiT) 
displays a more complex behavior with T as it depends on the backscattering amplitude, g\, which 
scales down under RG transformation and eventually behaves as g\ (T) ~ 1/ log T. We show, 
however, by direct perturbative calculations that Cs{T) is strictly linear in T to order gf as it 
contains the renormalized backscattering amplitude not on the scale of T, but at the cutoff scale 
set by the momentum dependence of the interaction around 2kF. The running amplitude gi{T) 
appears only at third order and gives rise to an extra T/\o^T term in Cs{T). This agrees with 
the results obtained by a variety of bosonization techniques. We also show how to obtain the same 
expansion in gi within the sine-Gordon model. 

PACS numbers: 



I. INTRODUCTION 

The hallmark of a Fermi liquid is the linear depen- 
dence of the specific heat C{T) on temperature. A devi- 
ation from linearity at the lowest temperatures generally 
implies a non-Fermi liquid behavior. This generic rule 
is satisfied in dimensions D > 1, e.g., the non- Fermi- 
liquid behavior near quantum critical points is character- 
ized by a divergent effective mass and sub-linear specific 
heat. On the other hand, the behavior of the best stud- 
ied non-Fermi liquids - one-dimensional (ID) systems of 
fermions - is more subtle. A ID system of fermions can 
be mapped onto a system of ID bosons. As long as these 
bosons are free, i.e., the system is in the universality 
class of a Luttinger liquid, the specific heat is linear in 
T despite that other properties of a system show a man- 
ifestly non-Fermi-liquid behavior. However, backscatter- 
ing and Umklapp scattering of original fermions give rise 
to interactions among bosons. If these interactions are 
marginally irrelevant, C(T) may acquire an additional 
logT dependence. 

In series of recent publications, several groups stud- 
ied specific heat of interacting Fermi systems in dimen- 
sions 1 < £> < 3 1,2, 3,4, 5, 6,7, 8, 9T0]. These systems 
are Fermi liquids, and the leading term is C{T) = 7T. 
The sublcading term is, however, non-analytic: it scales 
as AjjT^ (with an extra logT factor in 3D), and in 
1 < 13 < 3 the prefactor is expressed exactly via the 
spin and charge components of the fully renormalized 
backscattering amplitude^i^iiiS 
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where Ud is a number [02 — 3C(3)/(27r)], and /s(7r) and 
/c(7r) are components of the backscattering amplitude 
f{9 — -k) {9 is the angle between the incoming momenta). 
The spin and charge contributions to the specific heat 



can be extracted independently by measuring the specific 
heat at zero and a finite magnetic field (a strong enough 
magnetic field (J-bH ^ T reduces the spin contribution 
to 1/3 of its value in zero field). 

As D ^ 1, becomes T, and the universal sublcad- 
ing term in the specific heat becomes comparable to the 
leading term. In addition, the spin component of the 
backscattering amplitude in ID flows under a renormal- 
ization group (RG) transformation, and, for a repulsive 
interaction, which is the only case studied in this pa- 
per, scales as 1/logT in the limit T ^ [IH. The 
charge component, /c(7r), on the other hand, remains fi- 
nite. Judging from Eq. ([1]), one might then expect that 
the charge component of the specific heat in ID scales as 
T, while the spin component, Cs{T), scales as T/log^T. 

This simple argument is, however, inconsistent with 
recent result obtained by Aleiner and Efetov (AE) 
for the model of weakly interacting electrons. They 
developed a powerful "multidimensional bosonization" 
method in which fermions are integrated out and the ac- 
tion is expressed solely in terms of interacting, low-energy 
bosonic modes. In ID, AE showed that Cs{T) behaves 
as T/log^(r) for T ^ (in disagreement with the RG 
argument), and that the logarithmic flow of fs shows up 
in C(T) only at fourth order in the interaction. Simi- 
lar results have been previously obtained for the Kondo 
modeli^ and XXZ spin 1/2 chaini^, which are believed to 
be in the same universality class as ID fermions with re- 
pulsive interaction. [Earlier perturbative studies of C(T) 
in ID yielded different results: in Ref. [l3|, C{T) was 
argued to be linear in T to all orders in the interaction, 
whereas Ref. found that Cs(T) scales as T/log^T, 
both results are in disagreement with the result by AE.] 

The functional form of Cg (T) is not a purely academic 
issue. In a strong enough magnetic field ^bH ^ T the 
logT term is replaced by the logH one. Measuring the 
field dependence of C (T, H) , one can explicitly determine 
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the functional form of C{T,H = 0). We note in pass- 
ing that the issue of universal temperature corrections to 
thermodynamic quantities is not restricted to the specific 
heat. Number of researchers studied the universal tem- 
perature and wavevector dependence of the spin suscep- 
tibility — . Another example of a universal, non-analytic 
behavior is the T^/H behavior of the specific heat of a 
2D d— wave superconductor in a magnetic field^^. 

Absence of the logarithmic renormalization of Cs{T) 
below fourth order of perturbation is a rather non-trivial 
result in view of Eq. ([T]), but even more so because 
backscattering in ID contributes to the specific heat al- 
ready at first order in the interaction (see Sec lIII Bj l. 
Moreover, both first and second-order contributions to 
Cs{T) can be straightforwardly obtained in a computa- 
tional scheme in which they appear as contributions from 
low energies, of order T. In the RG spirit, one might ex- 
pect these terms to contain the running backscattering 
amplitude at a scale of order T. However, in ID, the exis- 
tence of a particular computational scheme, in which the 
answer comes from low energies, does not actually guar- 
antee that the corresponding coupling is a running one, 
as ID systems with a linear spectrum are well-known to 
exhibit anomalies, similar to Schwinger terms in current- 
current commutation relations. 

From computational viewpoint, the anomaly-type con- 
tribution to C(T) can be equally obtained either as a 
low-energy contribution, or as a contribution from high 
energies, of the order of the cutoff. In the latter case, the 
corresponding coupling is on the scale of the cutoff, rather 
than T. One then has to explicitly evaluate higher-order 
terms to verify whether the coupling is a bare one or a 
running one. 

This running vs. bare coupling dilemma was discussed 
actively in the earlier days of bosonizatio n^^i^^i^^ , and 
is related to a more general issue of how to treat prop- 
erly the high-energy cutoffs in theories with linear dis- 
persionsiS. 

Our interest in the ID problem is three-fold. First, 
we want to understand which of the backscattering cou- 
plings entering C{T) are the running ones and which 
are the bare ones. We argue below that anomaly-type 
terms should be treated as high-energy contributions, for 
which the couplings are at the cutoff scale. The running 
coupling appears in C(T) due to non-anomalous contri- 
butions, which can be uniquely identified as low-energy 
contributions. Second, we would like to check directly 
whether the g-ology model is a renormalizable theory or 
not, i.e., whether the dependence of the ultraviolet cut- 
offs can be incorporated into a finite (and small) num- 
ber of renormalized vertices. Third, we want to estab- 
lish parallels between the direct pcrturbative expansion 
in the backscattering amplitude in momentum space, and 
the real-space calculations within the sine-Gordon model. 
In particular, we want to understand how anomaly-type 
contributions appear in real-space calculations. This has 
not been considered in earlier works2ii^i2^i^ for which 
the main interest was a search for a contribution with 
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FIG. 1: A model interaction potential with two cutoffs; Ab 
near q = and Ab near q = 2kF- 

the running coupling. 



A. Model 

We consider an effective low-energy model of ID 
fermions with a linearized fermionic dispersion near 
±fci?, Efc = ±vp{k =p kp), and with a short-range four- 
fermion interaction U (q) . We set a fermionic momentum 
cutoff at a scale Af (generically, comparable to the lat- 
tice constant) , and assume that fermions with energies 
larger than vp-^h account for the renormalization of the 
bare interaction into an effective one, which acts between 
low-energy fermions, and depends not only on transferred 
momentum, but also on two incoming fermionic mo- 
menta. We then use the g-ology notations^^, and intro- 
duce three dimensionless vertex functions, gi, 52 and 54, 
which describe scattering processes along the Fermi sur- 
face with zero incoming and 2kp transferred momenta, 
zero incoming and zero transferred momenta, and 2kp 
incoming and zero transferred momenta, respectively. 
At first order in the interaction, gi = U{2kp)/{2Trvp), 
92 — g4 — U{0)/{2ttvp). The effective low-energy model 
only makes sense if the couplings gt vanish long before 
the scale of Af , otherwise the low-energy and high-energy 
sectors could not be separated. A way to enforce this 
constraint, which we will adopt, is to assume that the in- 
teractions gi are non-zero only for transferred momenta 
(either around zero or 2k p), which are smaller than Af. 
Accordingly, we introduce two "bosonic" cutoffs: Ab, set 
by the interaction with the momentum transfer near 2k p 
{gi vertex), and Ab, set by the interaction with a small 
momentum transfer (32 and 174 vertices), and request that 
both are smaller than Af . More precisely, we assume that 

Af-Ab,Af-Ab»T. (2) 
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The model interaction is shown in Fig.[T] We will see that 
there is an interesting dependence of the specific heat on 
the ratio Ab/Af, but no dependence on the ratio Ab/Af. 

The two-cutoff model with Ab < At has been used in 
the canonical ID bosonization approach, and in the sub- 
sequent analysis of the sine-Gordon model, ft was also 
considered in Refs|T^[l3 in the analysis of the electron- 
phonon interaction in ID. We note in passing that, to 
our knowledge, it has not been explicitly verified that 
the specific heat for the effective low-energy model is the 
same as for the original model of fermions with parabolic- 
type dispersion and a generic interaction U{q), i.e., that 
all contributions to C (T) from fermionic energies exceed- 
ing At can be absorbed into the three couplings gi. We 
also note that, in the bosonization procedure invented 
by AE, which is not based on the g— ology model, the 
cutoff imposed by the interaction is less restrictive that 
the fermionic cutoff (i.e., Ab ^ At), because in their the- 
ory the propagators of long-wavelength bosonic modes 
are obtained by integrating independently over fermionic 
momenta linked by the interaction. AE, however, only 
focused on the truly low-energy terms with the running 
coupling, which should not depend on the ratio Ab/At. 

Vertices 171 and 52 in the g -ology model are related 
to the spin and charge components of the backscattering 
amplitude 

as /c(7r) = 2.92 - gi, /s(7r) = -gi. Vertex 54 is re- 
lated to the forward scattering amplitude /(O) as 34 = 
/c(0) ~ /s(0)- For generality, we extend the model from 
the SU{2) symmetric to anisotropic case, i.e., assume 
that all three vertex functions gi (« = 1, 2, 4) have differ- 
ent values gi\\ and gi±, depending on whether the spins 
of the fermions in the initial state are parallel or oppo- 
site. For anisotropic case, the spin component of the 
backscattering amplitude splits into the longitudinal and 
transverse parts, and we have 
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The forward scattering vertex is invariant under RG 
renormalization, but the backscattering vertices gi and 
g2 flo w^^'^^ . Keeping only the processes with momen- 
tum transfers in narrow windows near either zero or 2k p 



(the windows are much smaller than the cutoffs Ab and 
|5arAjj), we have 
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where L = log Ep /E, E is the running energy, and /? 
functions depend on the couplings .gi||,_L ^-nd .92||,±^ In 
the one-loop approximation. 



/3i|| =-(.9?||+.9?±), Pi± 
I32\\ = -.9?|h (32± 



■9^. 



-25i_L (ffiii - .g2ll + 52_l) , 
(6) 



Re-expressing the couplings in terms of spin and charge 
components of the backscattering amplitude, we find 
that the spin amplitudes and fs± flow to zero un- 
der the RG transformation. The charge component of 
the backscattering amplitude /c(7r) = {g2\\ + 52_l) — .9i||, 
however, does not change under the RG flow. 

For the SU(2) symmetric case, /3i = —25^/(1 — 
91)1 P2 — /3i/2, and gi{L) renormalizes to zero as 1/L, 
while g2{L) tends to a constant value of a half of the 
charge amplitude, which is invariant under RG. 

As we said earlier, the key interest of our analysis is to 
understand at which order within the 5— ology model the 
running couplings appear in the specific heat, and what 
are the contributions to the specific heat which contain 
bare couplings. 



B. Results 

We first catalog our main results, and then present cal- 
culations in the bulk of the paper. We computed C{T) 
in a direct perturbation theory, expanding in powers of 
the couplings gi to order g^. To first two orders in gi, 
we found that the specific heat is expressed via bare cou- 
plings g2 and 174, and the effective backscattering cou- 
pling gi. At third order, we found an extra contribution 
to C(T), which comes from low-energies and contains 
a cube of the running backscattering amplitude on the 
scale of T. Explicitly, for the anisotropic case, we found 
for T <C Ab and neglecting 0{g^) contributions with non- 
running couplings 



C(T) = 



27rT 
3vF 
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1 + (5111 - 54|l) + (5111 - 54|l) + 54J_ + 2 ((-9211 + 92±) - 2.9 



?2||5l 



+ (.9?ll 



.9?J 



3.9?^(r)gi| 
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Here 54 and g2 are the bare couplings, and g^^ and gi± 
are the effective couplings on the scale Ab. The couplings 



gi\\{T) and gi±{T) are running couplings on the scale of 
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FIG. 2: One- loop diagrams for the interaction vertices. In the 
RG regime (external momenta are much smaller than A;,), the 
renormalizations of gix and giy are given by diagrams a), b), 
and d), while the renormalizations of 172^ and g2\\ are given 
by diagrams c) and d). All diagrams give rise to logT terms. 
For external momenta of order As, only diagram a) gives rise 
to the logarithmic term log(Af/Ab), as the two fermions in 
the particle-hole bubble can have momenta in the whole range 
between Ab and Af. For all other diagrams, the interaction 
constraints the internal momenta to be of the same order as 
the external momentum, and there is no momentum space for 
the logarithm. 



T - these are the solutions of the full RG equations, Eq. 
([5]), with effective giy and gi±^ serving as inputs. 

The effective couplings and gu^ differ from bare 
5i||,-L due to RPA-type renormalizations by ikp particle- 
hole bubbles made of fermions with momenta between 



Ab and Af. This renormalization comes from diagram 
a) in Fig. [5]). There are no such renormalizations for §2 
couplings, which retain their bare values. We obtain 



.9111 = .9111 - (.9f|| + 9n_ ) Lb+\9t\\ +^9i\\9n_ ) ^^(Sa) 
fiW + .9?± = .9' I + 5?± - 2 (g% + igngli) ib(8b) 

where Lb — log (Af/Ab). 

We emphasize that an RPA-type renormalization is not 
equivalent to RG, so that effective giy and ^ij^ differ from 
the solutions of (O, (O already at one- loop order. The 
difference is due to the fact that in the one-loop RG equa- 
tions, the RPA and ladder- type renormalizations of 171, 
and the ladder renormalizations of 32 are all coupled, 
while only the RPA diagrams lead to Lb terms in the 
renormalization from gi to gi . 

Coupling between the RPA and ladder renormaliza- 
tions in the RG regime is absent for the isotropic, SU (2) 
symmetric case. Then gi = gi/(l + 2giL\j) becomes 
equivalent to one-loop RG. Furthermore, in the symmet- 
ric case, the running gi{L) at the lowest energies behaves 
in the one-loop approximation as .9i/(l + 2giL). For the 
specific heat we then obtain 
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3.9? 
(1 + 2.giL) 



(9) 



The results of our direct perturbative analysis are in 
agreement with the results for the Kondo problemi^and 
XXZ spin chain^^ - for both models, the specific heat 
was shown to behave as T/ log'^ T at the lowest temper- 
atures. These two models are argued to be in the same 
universality class as the model of interacting electrons 
with the interaction in the spin sector. The same behav- 
ior was found by Cardy^ and Ludwig and Cardy^Z in 
their study of a conformally invariant theory perturbed 
by the marginal perturbation from the fixed point (the 
sine-Gordon model belongs to this class of theories), and 
by AE in their "multidimensional bosonization" analysis. 
In all these theories, the focus was on the universal terms 
which are confined to low energies i.e., are not anomalies. 
If only such terms are included, the full spin contribu- 
tion to the specific heat scales as T/ log'^ T in the SU (2) 



isotropic case, i.e., the spin part of the specific heat coef- 
ficient vanishes at T = 0. Our direct perturbation theory 
reproduces the same universal behavior in the spin sec- 
tor, but also generates extra contributions to the specific 
heat which contain effective interaction on the scale of 
Ab. 

To make the comparison with the bosonization and 
sine-Gordon model explicit, we re- write our result via 
spin and charge velocities vpUp and vpUa obtained by 
diagonalizing the gradient part of the Hamiltonian: 

^^p = (l + .9411 + .94± - .9111 (.9211 +-92-L - .9111)^, 

^^ = (l+.94ll - .94_L - gill (.9211 - .92_L - .9111)^(10) 

Using pO|) . one can re- write (O as 



C{T) = — ^ + ^ + _.g2 + 52 T gill T . 

I 



The last term in pT|) is the universal contribution from heats of two gases of free particles with the efi^ective ve- 
low energies. The first term is the sum of the specific 



locities Up and Wo-, which are the same as in (|10p except 
for (7]^ 1 1 and gi±^ are now the effective, renormahzed ver- 
tices. The term in the middle is an additional contribu- 
tion from the spin channel. Very likely, this contribution 
can be absorbed into the renormalization of spin velocity 
M(T — > U(T — ffij^, i.e., the specific heat can be re-expressed 
as the sum of the contribution with running couplings, 
and the specific heat of two ideal gases of fermions with 
bare charge velocity (albeit with gi), and the renormal- 
ized spin velocity. 

In the rest of the paper, we present the details of 
our calculations. In Sec. IIIB) we outline the compu- 
tational procedure, calculate the first-order diagram for 
the thermodynamic potential, and demonstrate explic- 
itly the sensitivity of the result for C{T) to the ratio 
of the cutoffs. In Sec lII CI and IIIDi we compute second 
and third order diagrams for the thermodynamic poten- 
tial, and discuss the fourth-order result. In Sec lIIIl we 
analyze the specific heat in the framework of the sine- 
Gordon model. Sec. IIVI presents the conclusions. Some 
technical details of the calculations are presented in the 
Appendices. 
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FIG. 3: First order diagrams for the thermodynamic poten- 
tial. Here and in the rest of the figures, the dashed line repre- 
sents the interaction. Diagram lb) does not contribute to the 
temperature dependence of the thermodynamic potential. 



S(0) =-TY, M = ^ + const, (14) 

such that C(o)(r) = 2t:T/3. 



B. First order diagrams 



II. PERTURBATION THEORY AND THE ROLE 
OF CUTOFFS 

A. Preliminaries 

In this and the next two Sections we set Ui? = 1. We 
restore vp in the final formulas for the specific heat. 

The specific heat of an interacting system of fermions 
can be extracted from the thermodynamic potential S 
via C{T) = -Td'^E/dT^. The thermodynamic potential 
is given by the Luttinger-Ward formula: 

E = S(«)-2T^ / I [log (GoG-) - EG + E i^E.G 

(12) ■ 



where 



77(0) 



E[ dk 
It. 



2 log (4 + 



(13) 



is the thermodynamic potential of the free Fermi gas per 
unit length Gq = {iuj„i ^ Cfe) ^, ffc is the dispersion, G = 
{iojm — e/c + S) ^, E is the exact (to all orders in the 
interaction) self-energy, and S^, is the skeleton self-energy 
of order ly. The skeleton and full self-energy are related 
via 5] = and are evaluated at finite T. Expanding 

both G and in Eq. ifTTj) in powers of the interaction, 
one generates a perturbative expansion for S in the series 
of closed diagrams with no external legs. 

The free-fermion expression for G(T) is obtained from 
Eq. I|13p . At low T, the momentum integration is con- 
fined to fc w ±kp and yields 



At first order, the T dependence of S comes from the 
bubble diagram crossed by the interaction line [diagram 
la) in Figl3]). This diagram contains two contributions: 
one with a small momentum transfer and another with a 
momentum transfer near 2k p- As spin is conserved along 
the bubble, the corresponding coupling constants are (74|| 
and respectively. 

The safe way to evaluate the diagram is to sum over fre- 
quencies first, as the frequency summation is constrained 
neither by the interaction nor by the fermionic band- 
width, and then integrate over the fermionic momentum 
k and bosonic, transferred momentum q. We will mea- 
sure g as a deviation from zero for (74|| term, and from 
2kp for (7i|| term, and, as we said, will cut interactions 
at \q\ = Ab for forward scattering process and at Ab for 
backscattering process. We linearize the fermionic dis- 
persion near the Fermi surface and set the cutoff of the 
integration over k at \k zL q/2\ < Af , Af > Ab, Ab. 

For small momentum transfer, we then obtain 



(1) 2g4|| f^^ ^£-9/2 cosh^ , ^ 

TT .In .In cosh 7^ + cosh ^ 



and for the momentum transfer near 2kp, we obtain 



Jo JO 



Ab /.Af-g/2 

dk- 



cosh y 



cosh ^ + cosh 



(16) 



Subtracting T-independent terms in (|15p and (|16p and 
introducing rescaled variables x = k/T, y ~ q/{2T), we 
re-write (fTSt and (fTHll as 



TT 



Jo Jo cosh y + cosh x 
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and 



and yields 



"(1) 



Ab/2T 

dy 



^ cosh X 
ax- 



coshi/ 



cosh X + cosh y 



(18) 



We immediately see that the first-order contribution to 
the thermodynamic potential vanishes if we formally ex- 
tend the integrals over x and y to infinity. Integrals (|17p 
and IjlSp are similar to the integrals which give rise to 
anomalies in the field theory22i. The integrands are odd 
under the interchange of x and y; therefore universal, 
cutoff independent contributions apparently vanish, but 
the 2D integrals are ultraviolet divergent, if we set T to 
zero. A finite T then sets an ultraviolet regularization of 
the divergent 2D integral and gives rise to finite terms 
in S which do not explicitly depend on the cutoff. By 
analogy with the field theory, hereafter we refer to these 
terms as "anomalies". 

For definiteness, we focus on the 2fci? contribution. 
Since At > Ab (in the sense of Eq.([2|), the integration 
over y extends to a much narrower range than that over 
X. In this situation, the most natural way to evaluate the 
thermal part of 'e!'^}! is to re-express ([T8|) as 



Ab/2T 



dy cosh y 



dx 



cosh X -I- cosh y 



(19) 



The integral over x now converges, and, because Af > Ab, 
we can safely set the upper limit of the x integral to 
infinity. The x integration then can be performed exactly 



"(1) 



A /•Ab/2T 

- .9111^^ / 

71" Jo 



dyy cothy, 



.9111 ,2 2 



2tt 



Ab/2T 

I dyy{cothy-l). (20) 





The thermal part of 'e!'^'^ 



comes from the second term 



41 = const -^31 II 



(21) 



Observe that the piece is independent of the cutoff. 
Furthermore, in this computational procedure, the fre- 
quency sums and the momentum integrals are fully ul- 
traviolet convergent, and Eq. pT|) comes from small mo- 
menta fc, g ~ T <C Af , Ab. 

Alternatively, however, we can evaluate the integrals 
in (|18p by integrating over the (dimensionless) bosonic 
momentum y first. To do this, we neglect y in the upper 
limit of the integral over x (we will check a'posteriori that 
this is justified), and re-express as 



4.91 



Ai/T 



dx cosh X 



Ab/2T 



dy 



(22) 

/ Q cosh X + cosh y 
It is tempting to set the upper limit of the y integral to 
infinity, as this integral converges. However, one has to 
be cautious as there is a range of x where cosh x > cosh y 
for any y. To see how this affects the result, we represent 
the y integral as 



J 



(■Ab/2T 



dy 



dy 



dy 



1 



dy 



Jo cosh X -I- cosh y Jo cosh x -|- cosh y jAb/2T coshx + cosh y sinhx coshx J(A|^/2T) 1 + e** ^ 

I 



(23) 



We replaced cosh x and cosh y in the last term by the ex- 
ponentials, as y are large, and we anticipate typical x to 



be large as well. The remaining integration is straight- 
forward, and we obtain 



J 



"(1) 



const -f -gill - / dx log f 1 + p/^'^''/^^ 

3' " TT Jo ^ 



(24) 



The first term is the contribution from low energies ~ the 
same as in ((2T|) . but with the opposite sign. The second 
term by construction is the contribution from energies 
much larger than T. Evaluating the second integral, we 



find that it also contributes a term to ^^2kF ■ 



4gi 



'^-^T^ j^''^ dx\og(l + e^~^-'^'^) 



27r 2 
const + —gi\\T 



(25) 
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The term in ((251) comes from x ^ Ab/2T. It is es- 
sential that these x are smaller than the upper limit of 
x— integration, otherwise such a contribution would not 
exist. Substituting this back into we find that the 
high-energy term is opposite in sign and twice larger than 
the low-energy one, so that the sum of the two contribu- 
tions is given precisely by Eq. l|2ip . Going back through 
the derivation of (j25p . we see that typical y and x are near 
Ab/2T, well below the upper limit of the x integration. 
In this situation, the neglect of y in the upper limit of 
the integral over x is legitimate, to accuracy exp — Ab/T. 

We see therefore that ^2^^^ can be equally well obtained 
either as a low-energy contribution or as a high-energy 
one. This is a hallmark of an anomaly. The same is true 
also for the forward scattering term S^^q: the term 
again can be equally obtained as a low-energy contribu- 
tion or as a contributions from energies of order Ab- 

Combining the results for backscattering and forward 
scattering, we obtain for the specific heat 

= 1^(3111-3411). (26) 

As we said in the Introduction, the gi\\ term in (|26p 
is not present in the standard bosonization approachii. 
The physical argument is that the gi\f^ should only ap- 
pear in C(T) in the combination gi\\ ~92\\ as the two ver- 
tices transform into each other by interchanging external 
momenta without interchanging spins, and therefore are 
physically indistinguishable^^'"^". The first order diagram 
with g2\\ is a Hartree diagram with two bubbles connected 
by the interaction at exactly zero transferred momentum 
[diagram lb) in Fig[3]. As each of these two bubbles rep- 



resents a total electron density, this diagram obviously 
does not depend on T. By the argument above, the dia- 
gram with 3111 also should not depend on T. This consid- 
eration is, however, only valid if the cutoffs are infinite. 
For finite cutoffs, there appears an extra "anomaly-type" 
contribution, in which 31 1| appears in the combination 
with 32||, as we just demonstrated . A similar reason- 
ing within real-space consideration has been presented 
in^"^. Another argument for presence of the 31 1| term is 
based on the observation that fermions with the same 
spin do not interact via a contact interaction; hence the 
interaction should drop out of the results in this limit. 
This implies that the observables, such as C(T) must 
depend separately on the combinations 34] | — 3i|| and 
g2\\—gi\\ — Eq. (pS)) is consistent with this argument 
as the limit of a contact interaction, i.e, for 3]^|| = 34||, 
C(i)(T) vanishes. 

The interplay between low-energy and high-energy 
contributions to 2 can be also understood if one inter- 
changes one momentum integration and one frequency 
summation and expresses 'E,^^^ via the polarization bub- 
ble as 

4Z = -9i\\Tj2 r dqli2kA<l,^)- (27) 

The sub-indices indicate that the momentum integration 
is confined to q near zero or near 2kp. 

For briefness, we consider only the backscattering 
term. The polarization bubble 112^^ (g, il) is given by 



Il2kAl.^) = 2TY, j ^GR{uj + n,k + q)GL{uJ,k)+GL{io + n,k + q)GR{Lo,k) (28) 



27r r"^~4A2 V('?-*^^)^-4fc2 ' {q + iVlf-Ak 



T^zi^og-^-sl dkknF{k){ ^^ ... + ..J ). (29) 



where GL,R{kuj) is the Green's function of right/left mov- 
ing fermions and npix) is the Fermi function. The first 
term in Eq. (|29p is the zero-temperature Kohn anomaly, 
the rest is the thermal contribution. The integration over 
k gives the result for Il{q, il) in terms of di-Gamma func- 
tions'^^, but for our purposes it is more convenient to use 

Substituting ^ into we find 

4Z^~9i\\{Q + P), (30) 

where Q and P are the contributions from the Kohn 
anomaly and from the thermal piece in (j29p , respectively. 



The temperature -dependent part of the Q term is 

comes from low energies regardless of whether the sum 
or the integral is done first. In both cases, we get, up to 
a constant, 

Q = Ql^ -^TV3, (32) 

where subindex L specifies that this is a contribution 
from low energies: fl,q ^ T. The second, thermal, term 
is determined either by low or by high energies, depend- 
ing on the order. If the momentum integration is done 
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first, the non-zero result is obtained only because Ab is 
finite; otherwise, the integration contour can be closed 



in that half-plane where the integrand has no poles. Re- 
arranging the integrals, we rewrite this contribution as 



J 



P 



^ /■Ab roc f Y Y 



g ^ roc roc 

—Ty / dq dkknp{k) 
^ ^Ja, Jo 



1 



1 



{q - - ik^ (q + in)^ -4k'^J ' 
I 



(33) 



As the Fermi function in confines the fermionic mo- 
mentum to fc ^ T, and q > Ab is large, we can neglect 
4fc^ compared to (q ± ifl)'^ in the denominator. This 
simplifies P to 



dq 



Ab 



1 



1 



{q-iny {q + iny^ 



Performing the momentum integration, we obtain 

^2 — Ab 



P 



4:TTT 



A2 



(34) 



(35) 



Evaluating the integral, we find that P does not depend 
on the cutoff, and equals to 



P = Ph = 



2ttT' 



(36) 



where subindex H specifies that this is a contribution 
from high energies fl,q ^ A^,. 

Alternatively, P can be evaluated by doing frequency 
summation first. One can easily check that, to order T^, 
the frequency sum can be replaced by the integral. The 
frequency integral is non-zero only for q < 2k, otherwise 
the poles in are located in the same half-plane, and 
the frequency integral vanishes. Evaluating the frequency 
integral and then the integral over q, we reduce P to 



P^Pl 



dkknpik) 



(37) 



where subindex L specifies that this a contribution from 
low energies VL ^T. We see that Ph = Pl, i-e-, the same 
result for P can be obtained either as a high-energy con- 
tribution, or as low-energy one. In both cases P is for- 
mally independent of the cutoff, and the total backscat- 
tering part of S^^^ is given by 



^2fcs- = const — [Q + P) = const - 



(38) 



which coincides with pT|) . 

Which of the two ways (low-energy or high-energy) is 
physically correct? As we discussed in the Introduction, 
if 'e!^!^ comes from low energies (of order T) , one should 



expect T log T terms in C (T) already at the next (second) 
order; on the contrary, if it comes from high energies, no 
such terms are expected. AE suggested implicitly that 
the correct procedure is to take the average of two possi- 
ble orderings, i.e., to represent frequency summation and 
momentum integration in (|27p as 



(39) 



In this procedure, P in Eq. 

1 



P 



tPl 



33)1 is a sum 
1 



H, 



(40) 



with Ph ^ Pl = Total S^^^^ is the sum of P 

and Q [see Eq. (|38|) ]. where Q = Ql = — 7rT^/3 comes 
from low energies. Adding P and Q, we find that the 
low-energy contributions cancel out, and the net result 
for 'B^^kp high-energy contribution. 

Eq. ()39l) . however, contains some ambiguity, as one 
could equally well can re- write P as uPl -\-{\ — a)PH with 
an arbitrary coefficient a. Then, the balance between the 
low and high energy contributions to S*^^-* would depend 
on a. Whether S^^^ contains running or bare coupling 
gi\\ can only be established by an explicit computation 
to the next (second) order. This is what we will do in 
next Section. 



C. Second order diagrams 

The second-order diagrams for the thermodynamic po- 
tential are shown in Fig. 31 There are two different types 
of diagrams, obtained by inserting either self-energy cor- 
rections or vertex corrections into the first-order diagram. 
The diagram with self-energy insertions [diagram 2a)] is 
readily computed either explicitly, or by evaluating first- 
order self-energy and substituting the result into the first- 
order diagram. We will not discuss computational steps 
(they are not qualitatively different from those to first 
order) and present only the final result: the diagram 2a) 
yields a regular contribution to S of the form 



7(2) 
^2a 



37r 



(5111 



5411 



(41) 
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FIG. 4: Second order diagrams for the thermodynamic po- 
tential. 

Next, there are vertex correction diagrams 2b) and 2c), 
which involve the forward scattering vertex 174 (small 
transferred momentum and 2/ci? total incoming momen- 
tum) and §2 vertex (small transferred and small total 



incoming momentum). The contributions of order g^p 

94±^ 5211' 52_L' order .g2||5i||j a-rc all expressed 

via the bilinear combinations of the bubbles for right and 
left movers 

It: q 

whose sum is the total polarization bubble 
Iiq=Q[q,^) — IlflXg, il) + ni((7, il). The evaluation 
of ^En / d<l^Mi ^^)nj(g, ^) = L, R) is straightfor- 
ward, and the result does not depend on the order of 
momentum and frequency integrations. We have, up to 
T-independent terms. 



3 ' 



/ dq[lll{q,n)+nl{q,n)] =-^tY^ f d<Z^^5^ = -^T^|r!|-const 
I dqIl,{q,n)nn{q,n)^-^TY^ I ^9fi2^ (43) 



These expressions give rise only to regular terms in the 
thermodynamic potential and, consequently, to T terms 
in the specific heat. Collecting combinatorial factors, we 
find that the contribution from cancels out among 
diagrams 2b) and 2c), while the rest yields 



7(2) _ 



~3^ 



2 

.94± 



.92,1 



2 

32± 



3211.9111 



(44) 



Non-trivial second-order contributions are associated 
with the vertex corrections due to backscattering ampli- 
tude in diagram 2b). These are most easily expressed via 
the square of the 2kp polarization bubble as 



^2kf 



tE/^9(.9?|| 



gi^U'iq,n). (45) 



The evaluation of T J dqllli,^ (g, fi), presented in Ap- 
pendix A, gives 

Tj2j dqlll,^ (g, fl) ^ ^ (1 - 2Lb) , (46) 

where, we remind, Ly, = log(Af/Ab). Combining (I45|) 
and we obtain 



^(2) 
^2kp 



6 



9h 



gi^ (1 - 2Lb 



(47) 



The logarithmic term in l|47p is the contribution from 
particle-hole bubble with fermions with energies between 



Ab and Af. One can easily verify that it coincides with 
one-loop renormalization of the backscattering amplitude 
(7i|[, which also comes from the bubble diagram. Indeed, 
according to ([5a|) . the effective coupling gi\\ on the scale 
Ab is, to order g^. 



5111 



.9111 



(5iii +5i_L)-^b- 



(48) 



This is precisely what one obtains by combining ^ 



(1) 
2kF 



from ([55]) and the logarithmic term in "B/^f!^ . We see that 
at low T, the effective giy is the renormalized coupling on 
the scale of Ab rather than on the scale of T. This implies 
that the correct way to interpret the anomaly in ^2^'^ is 
to treat it as a purely high-energy contribution. This 
agrees with the "symmetrized" procedure of Eq. ([55]) . 

If the fermionic and bosonic cutoffs differ significantly, 
i.e., Af 3> Ab, then gi\\ flows logarithmically in the en- 
ergy interval Ab <C <C Af. This flow freezes, however, 
at ~ Ab . We also emphasize that the non-logarithmic 
term in ^2kF independent of the ratio of the fermionic 
and bosonic cutoff, just like first-order gi\\ — g^\ term. 
Another way to see this is to adopt a different computa- 
tional procedure for the backscattering part of diagram 
2b. Namely, by virtue of 2kp scattering, two pairs of 
fermions from different bubbles have nearly equal mo- 
menta. Combining these two pairs into two bubbles with 
small momentum transfers and integrating independently 
over the two running momenta in these two new bubbles. 
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one can re-express the contribution via the product of 
two polarization bubbles with small momentum trans- 
fers. This procedure was employed in earlier work for 
D > 1^, and by AE for ID. It is justified, however, only 
when the momentum dependence of the interaction is 



J 



weak up to the fermionic cutoff, i.e., in a formal limit 
when Ab S> Af (which is opposite to what we assume 

here). Applying this procedure, one can re-express ^2kF 
as 



E^'\2kF) = -27r (.g?|| +.9?i)T^ f dqUL{q,mR{qM) ^ { 



all 



(49) 



This agrees with Eq. (|47)) without the logarithmic term. 

Combing the contributions to S^^' from Eqs. ()41l44p 
and (fiS)) . we obtain for the second-order specific heat 



J 



C(2)(T) 



2tiT 



(glw + gli) + ^ ((gill - gA\\f + .gL) + +^ ((3211 - .9111)' + +9li_ + qIa) 



ttT 



3vf 



3vp 



(50) 



Note that the bilinear combination of 32 and gi in the 
last term of (150]) is precisely the same combination of the 
backscattering amplitudes as for the non-analytic term 
in C(T) in higher dimensions (c.f. Eq. (1)): 



(5211 - ffiii)^ +g2_L + 



in) 



(51) 



We see that besides the logarithmic term which trans- 
forms .gi|| into the specific heat C'^^(T) also con- 
tains the "universal" second-order terms that do not de- 
pend on the cutoffs. This poses the same question as 
before - are those couplings the running ones (on the 
scale of T) or the bare ones (on the scale of a cutoff)? 
On one hand, the combination of the second-order g2 
and gi terms in C'^'^^T) is the sum of the squares of 
charge and spin components of the backscattering am- 
plitude, Eq. ((5T|) . As the spin amphtude flows un- 
der RG and acquires logT corrections, one could expect 
TlogT terms at the next, third order. On the other 
hand, all constant terms in C'^^^T) can be formally 
represented as non-logarithmic renormalizations of gi|| 
and 3411 . This renormalization involves the static bubble 
Uq^oiq,^ = 0) = TX;^/ dkG{k,uj)G{k + q,uj), which 
is an anomaly by itself-it can be viewed as coming from 
low-energies (of order g), if we sum over uj first, or from 
high energies, of order Af, if we integrate over k first. 
It is then unclear a priori whether the scattering ampli- 
tudes in C(2)(T) are the amplitudes on the scale of order 
T or on the scale of the cutoff. To verify this, we need to 
compute explicitly third-order diagrams. 





3a) 



FIG. 5; Third-order diagrams that potentially give logarith- 
mic contributions. In diagram 3a), momenta k and p are 
counted from iikp, respectively. 



D. Third order diagrams and beyond 

1. Third order diagrams 

We analyze the third-order diagrams in two steps. At 
the first step, we analyze possible logarithmic terms in 
S at third order, searching for logT terms, and also for 
terms which contain logAf/Ab. We will show that there 
are no logT terms at third order, whereas all logAf/Ab 
can be accounted for by renormalizations of the <7i||,± 
vertices. At the second step, we will show that there 
exists a universal third-order term which starts to fiow 
at the next, fourth order. 

a. Logarithmic contributions. Diagrams that poten- 
tially contain logarithmic contributions are shown in 
Fig. [51 Diagram 3a) is expressed via the cube of the 
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polarization bubble at 2fci?: 



for diagram 3b) is 



"(3) 
^3a 



(52) 

The computation of T J2n J ^l^lkp ^) lengthy, and 
we present it in Appendix B. The result is 



TJ2 [ dqIll,Jq,n)^^{Ll~U) 



(53) 



"(3) 
"36 



2nT' 



m\-m 



9l\\ 



Lb 



(55) 



Diagram 3c) in Fig. [5] is a vertex renormalization of 
the second-order diagram 2c) in Fig. 01 One of the lines 
of the second-order diagram is gi and the other one is 92 ■ 
Inserting the 2fci? bubble into the gi line, we obtain for 
diagram 3c) 



where, we remind. Lb = log (Af/Ab). AU potential log^ T 
and logT terms cancel out, and the only logarithmic de- 
pendence left involves the ratio of the cutoffs. Substitut- 
ing ([53|) into ([52| . we obtain 



77(3) 

'—'3a 



(.9?|| +39i\\9l±) {lI - Lb.) 



(54) 



Diagram 3b) is a vertex renormalization of the second- 
order diagram 2a) in Fig. 31 The vertices in diagram 
2a) can be both 171 or one of them can be gi and the 
other one g^. The ikp bubble in diagram 3b) of Fig. [5] is 
inserted into the gi line in both cases. The total result 



"(3) 
"3c 



-5211 



2 

9i± 



(56) 



Note that only 2k p couplings in C(T) are renormal- 
ized. The g2 coupling in C(T) remains at its bare value. 
A renormalization of g2 could potentially come from dia- 
gram 3d) , but this diagram contains no Lb terms because 
all internal momenta in this diagram cannot deviate from 
external momenta more than by Ab, i.e., there is no space 
for the logarithm in momentum integrals'^^ . Therefore, 
the logarithmic part of the third-order specific heat is 
obtained by combining the results from Eqs. (|54l55p and 



9l, 



35i||.9i_L 



2ttT 



Lb 



9i\\ +3gi||.gi_L + 2 Ui|| -3411 



^9211 



)( 



9l\\ 



9l^ 



(57) 



r 



The third-order specific heat can be obtained from the 
results of the first and second orders by replacing the 
bare couplings by their renormalized values, 

In particular, the h\ term in Eq. ([Ff|) accounts for the 
third-order ladder renormalization of g\\\ {g\\\ {g\^\^ -\- 

3g^_i_gi||)L^, see (|8a|) ] in the first order specific heat 
[Eq. ([26]) ]. The Lb terms account for the renormaliza- 
tions of the g^w-, and g\^ -I- g\^ terms in C^^^(r) 
[Eq. ipS]) ] according to 



ffill 
2 

ffill 

2 I 2 

fill + .9l_L 



[see Eqs. (|5al8b)) ]. 



+ .9?_L ) 



-2gi|| (ffiii +g?_L) Lb 



-2 .9^11+ 3.91 II. 9^^ Lb, 



h. Universal contributions. To this end, we have not 
obtained a term with the running coupling on the scale 
of T. We now demonstrate how such a term is generated 
at third order. To do this, we compute the constant, 
cutoff-independent term in S'^' . We will not attempt to 
calculate this term using Eq. ([5^ . as the calculations are 
quite involved. Rather, we assume, by analogy with the 
second-order calculation, that this constant term is inde- 
pendent of the ratio of the cutoffs and can be evaluated 
in the same computational procedure as the one that led 
us to Eq. (|49|) . i.e., by reducing the 2kF problem to the 
small 9 one, and representing the third-order diagrams 
as the products of two triads. The same procedure was 
employed by AE. 

The relevant diagrams here are diagrams 3a) and 3d). 
Using the triad method, we obtain for their sum 



^(3) ^^(3) 
'-^sum "3a 



"(3) 
"3(i 



-.911 



/ / rf92n3(9i,92,^^i,f^2)n3(-9i,-92,rJi,rj2), 

ill 02 



(58) 
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where 



113(91, 52, f^i,i^2) / dkGR{k,uJk)GR{k + qi,ujk + i^i) GR{k + q2,uJk + ^2), 



(59) 



The integration in ([59]) is straightforward, as aU integrals 
converge, and we have 



113(91,92,^^1,^^2) 
1 



1 / i^l2 + 92 i^i + qi 



2tt \ifl2 — 92 i^i — 9i 



X 



+ ^2) - {qi + 92) 



(60) 



Substituting into 



we obtain 



-(3) ^ 2 

^sum = -5lll.9lJ. 



J dq, J dq2Tj2Tj2 



1 



^2 



1^2 -q2 iili - 91/ V*^2 + 92 ifli+qi 
1 



+ ^2) - (91 + 92) ii^^i + ^2) + (91 + 92) ■ 

I 



(61) 



The computation of the double momentum integral and 
frequency sum requires special care. The most straight- 
forward way is to sum over frequencies first, as the fre- 
quency sums are not restricted by cutoffs. Performing 
the summation, and using the symmetry between qi and 
92, we find after some algebra 



2£L = -9i\\9iJl / dqi / ^92 T 



(62) 



This obviously implies that the momentum integral is 
confined to high energies, of order Af, and does not 
contain a term. 

However, this is not the whole story. The new un- 
derstanding is obtained if we perform computations in 
different order, by integrating over momentum first. The 
computation is again lengthy, but straightforward, and 
yields 



77(3) 

"""sum 



x(l 



6nuoSn2fi)F{^i,^2,Af) 



(63) 



where 6a,b is the Kronecker symbol, and F(Oi,r22,Af) 
approaches a constant (= 1) when frequencies are much 
smaller than the fermionic cutoff Af . At frequencies com- 
parable and larger than the cutoff, F is rather complex, 
but the part of F relevant for our purposes is 



F(r!i,l]2,Af)-l-- 

TT 



Af 



A? 



1^2 



Af 



nj + Af 



(64) 



There are other terms in F, but they do not lead to a 

term in SsSL. 

The double frequency sum in (|63p then reduces to 



6 



rill 



Af 



(65) 



where the summation is now over all Matsubara frequen- 
cies, including fli,n2 = (we used the symmetry be- 
tween ill and ^12)- The sum T'^q^ ^Sfia ^ confined 
to large frequencies, and does not lead to term in 
si^-*. If Af were infinite, — would be the only out- 
come of (|65|) . For finite Af, one has to be careful as the 
second term in cannot be neglected for 0.2 > Af . Re- 
placing the sum over f22 by the integral, we obtain that 
the contribution from the second term in (|65p reduces to 
-(3/7r)rX]oi l^il = const + T^. Adding this result and 
the -T^ term in we find that the term in si^^ 
vanishes. This agrees with (|62p . However, we now see 
that the vanishing of the term in E^'' is the result 
of a cancellation between two physically different con- 
tributions. The —T^ term in (j65p is a truly low-energy 
contribution, which survives even if we set Af = 00. This 

term comes from fli — Q2 — 0, and from vanishingly 
small 9i, 92 in the momentum integrand. A very similar 
term leads to a non-analyticity in the spin susceptibil- 
ity ^. The coupling for this term, gf±gi\\, is then at the 
low-energy scale (~ T), and should be fully renormal- 
ized within RG. On the other hand, the compensating 

term comes from large energies, of order Af, and is 
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therefore a high-energy contribution. The corresponding 
couphng is then at the high energy scale, and it should 
remain constant under the RG transformation. 
As a result, sium becomes 

^ill = -^T^ {glALT)9i\\{LT) - ghgiw) ■ (66) 

where 51 (Lt) is the 2kF coupling on the scale of T, and 
gi without argument is the coupling at the cutoff scale. 
The low-energy contribution in (|66p coincides with the 
result obtained by AE (modulo a factor of 2). AE did 
not evaluate the high-energy contribution in (I66p . We 

did not attempt to obtain 2a for an arbitrary ratio of 
Ab and Af. We expect that the low-energy contribution 
is independent of the ratio of the cutoffs. At the same 
time, the high-energy term in (|66p may depend on the 
ratio of the cutoffs, i.e., the (—1) factor between low- 
energy and high-energy contributions in (j66p may only 
hold for Ab > Af, when the "triad" calculation is valid. 
In any event, the high-energy term in (j66p is a regular 
term and is therefore of little interest. 
For completeness, we also note that there exists an- 
other high-energy contribution of order gij_gi\\, obtained 
by inserting the first-order renormalization of the Fermi 
velocity into the second-order backscattering diagram. 
This contribution can be easily evaluated in the same 
way as (|66)) and yields 

^ixtra = -T^T^gl^dlW (67) 

If is independent on the ratio of the cutoffs, the 
high-energy terms in ([55]) and (p7)) cancel each other, 
i.e., the net result is only low-energy contribution. This 
cancellation is likely accidental, however. 

Assembling logarithmic and universal constant term at 
third order, evaluating the specific heat, combining with 
first and second-order diagrams, and using the RG flow 
of the couplings, we obtain the full result for the specific 
heat C(T), Eqs. Q and (H]). 

2. Fourth order diagrams 

For completeness, we also computed explicitly the 
fourth-order, four-bubble backscattering diagram for the 
thermodynamic potential We indeed found a log T 
term obtained by combining the "zero-energy" contribu- 
tion to Ssum [Eq. (|66p ] with an additional polarization 
bubble U2kp{0,0) esc logT. This one accounts for the 
logT renormalization of the running couplings gi±{LT) 
and (7i||(Lt) in We searched for possible other 

T^logT contributions, using the same method as at the 
end of the previous section. Namely, we assumed that 

log T terms must be independent of the cutoff ra- 
tio, set Ab > Af, and created two "quaternions" by as- 
sembling four fermionic propagators with close momenta 



k « kp, k + qi, k + q2, k + qs, and p « —kp, P + Qi, 
p + q2, P + \q\i ^ [see Fig [5]. We integrated inde- 
pendently over k and p in infinite limits (we recall that 




FIG. 6: Fourth-order diagram with four bubbles. 



this is possible only if the cutoff imposed by the interac- 
tion is irrelevant), then integrated over qi and summed 
over corresponding frequencies. We found logT terms 
from particular regions of frequency summations and in- 
tegrations over three bosonic momenta qi, (72, and q^; 
however, all such T^logT terms cancel out. Therefore, 
the only non- vanishing log T contribution at fourth 
order is the "zero-energy" one. All terms in S up to 
fourth order are anomalies, and corresponding couplings 
are at energies 0(A). 



III. COMPARISON TO THE SINE-GORDON 
MODEL 

A. Model 

A well-established way to treat the system of ID 
fermions is bosonization, which allows one to map the 
original problem onto the quantum sine-Gordon model. 
It is instructive to see how the results of the previous Sec- 
tions can be obtained within this model. We bosonize the 
operators of right- and left-moving fermions, Ra and Lq., 
in a standard way 

Ra{x),La{x) = i — exp [±i {(j)a{x) T 0a{x))] , a =t,i, 

where a is a short-distance cutoff related to the mo- 
mentum cutoff introduced in the previous Sections via 
a = Aj^. Upon bosonization, the part of the fermionic 
Hamiltonian parameterized by couplings g^ and 172 is 
mapped onto the Gaussian part of the bosonic Hamil- 
tonian, Hg = Hq'' -f Hq^ , where 
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H. 



^ j dx{l+ 3411 ± 5r4_L + ff2|| ± 52_l) {dx(t>p,a f + (l + 54|| ± 54± - 52|| T 52±) {dxd p,af , 



(68) 



r 



and the charge and spin bosons are defined as <j)p.a — 

(01 ± (j){) 1^/2 and 6p,a = (fi*! ± ^i) /v^- The 2fci. scat- 
tering, however, leads to non-hnear, cosine terms in the 
bosonic Hamiltonian. For a local, delta- function interac- 
tion, the cosine term only comes from 2fcF-scattering of 
fermions with opposite spins (coupling gn_). However, 



J 



for an arbitrary, non-local interaction, there is also a co- 
sine term which comes from 2fc_F-scattering of fermions 
with parallel spins (coupling gi\\)- Introducing a finite- 
range interaction V12 = V{xi — ^2), we map the 2fci? part 
of the fermionic Hamiltonian onto 



H 



1||,-L 



(2™)^ 



dxidx2Vi2 cos 



27r {4>p{xi) - (j)p {X2)) + 2kF [xi - X2) 
I 



/2^(0,(a;i)T 0.(3^2)) ■ (69) 



For a local interaction, V12 = Vii5{xi — X2), Eq. 
reduces to the usual sine- Gordon model. 

The universal term in the thermodynamic poten- 
tial (the analog of the universal term in Eq. ([66]) for 
the SU(2) symmetric case) was obtained by Cardy2^ and 
Ludwing and Cardy^ for a general case of a conformal 
theory perturbed about a fixed point by a marginally ir- 
relevant operator, and we just refer the reader to that 
work. The first and second-order terms in gi, however, 
have not been obtained explicitly in the sine-Gordon 
model before. Our goal is to demonstrate how the anoma- 
lous terms of order 171 and appear in the thermo- 
dynamic potential, and, in particular, how the gi cou- 
pling gets a logarithmic renormalization on a scale of the 
bosonic cutoff in this model. We will see that to get 
this renormalization, and also to obtain g^ term with a 
correct prefactor, one must consider a finite-range inter- 
action and keep the range of the interaction larger than 
the short-distance cutoff of the theory. 

The thermodynamic potential per unit length is given 

by 



-^log j D^e^p{~ [^G + ^i|| 



(70) 



where Sa, with a = G, 1||, and 1 _L, are the actions cor- 
responding to the Gaussian and 2kp parts of the bosonic 
Hamiltonian, respectively. Expansion in 5*111 -fS'ij^ gener- 
ates perturbation series for S. In the absence of backscat- 
tering (51 = 0), the bosons are free and theory is exactly 
solvable for arbitrary 52 and 54. One can then construct 
the perturbation theory in gin and gi± about the free- 
boson point. To make a connection with the previous 
Sections, however, we will perform the perturbative ex- 
pansion in all coupling constants rather than only in gi . 
This means that the averages generated by an expansion 
in 5*111 and 5ij^ will be taken over a free Gaussian action, 
Eq. dSHD with 34 = 52 = 0. 



B. First order 

The first-order term is obtained by expanding the ex- 
ponential in (|70p to first order in 5i||. Performing the 
averaging, we obtain 

S^^ll) = 2 / dxV {x) Ap{x, 0)Aa{x, 0) cos(2fci.a;), 

J\x\>a 



where 



(71) 



(72) 



As the averages are calculated over a free Gaussian ac- 
tion, Ap and Aa- are equal to each other and given byii 

1/2 

Ap{x,T) = A^{x,t) 



4 sinh TTxT + sin tttT 

The 2kp polarization bubble in the x, t space is 

1 

^2kp {x,t) = -Al ,^ {x,t) 



(73) 



4 sinh TTxT + sin tttT 
(74) 



Therefore, the first-order result reduces to 

S^^ll) = -2 y" dxV (x) H2fc^ {x, 0) cos{2kFx). (75) 

Note that this is nothing more than the first-order dia- 
gram la) in Fig. O written in the x, r space. Expanding 
H2fej^ {x,0) for X <C T^^, we obtain 



H2fe^ (a;,0) 



y2 



471^3 



4 sinh nxT 

The universal, constant term in Il2kp {x, 0) gives a T de- 
pendent part of Sj]^^^ 



^(i||) 



const- -5i II T^, 



(77) 
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where = (l/27r) / dxV (x) cos[2kpx]. We see that the 
bosonization result, Eq. ([77)) . agrees with the diagram- 
matic one, Eq. ( [2T|l . obtained for Ab < Af. This last 
condition is implicit in bosonization, as Eq. (|74p is valid 
only if the fcrmionic cutoff exceeds the bosonic one. No- 
tice that the final result Eq. ([77|) is formally valid also 
for a local interaction. However, the limit of a local in- 
teraction cannot be taken at the very beginning. Indeed, 
in this limit reduces to a constant and does not con- 
tribute to the T dependence of S. 

1. Gaussian form of Hi^^ 

The (7i|| term in the specific heat can be also obtained 
by reducing H^^ in Eq. (|69p to the Gaussian form, similar 



to what was done in Ref. \2§\ for spinless fermions. For 
completeness, we repeat this derivation here for fermions 
with spin. Indeed, can be written as the convolution 
of the 2kp components of the density 




(78) 

where p2k^Ax) = Ra (x) Ll (x) e^'''^^ + H.c = 
|^g2iV^F0„(x)g2i/cFX _^ jj.c.) /27ra. Performing normal or- 
dering in the product of two exponentials of the bosonic 
field and Taylor expanding the difference (j^a (xi) — 
4>a (^2) under the normal-ordering sign, one obtains 



g2iV^0a(a;i)g-2jV^()i„(2:i) g2i\/7F[0a (a;i)-0a (3:2 )] . g-j^-p 



1 - 27r {dx4>af [Xl - X2f' + 



-'i7,{<t>{c{xt)-4>^(x2)f) 



[xi - X2Y 

I — 



— c- number — 2Tia^ {dx(t>a) 



At the level of operators, Hiu then reduces to 



H 



111 



-5i|| 



Combining this result with the Gaussian part of the 
Hamiltonian, we obtain a new effective Hamiltonian 



The specific heat, corresponding to Hq, is given by the 
first term in Eg. pT]) . Notice that, in contrast to conven- 
tional bosonization which treats the gi\\ interaction only 
as an exchange process to the g2\\ interaction and, there- 
fore, contains only a combination of g2\\ — Kp „ in 
([75)1 contain two combinations: 17211 ~ 9i\\ and g^\ — gi\\. 



m 



where the charge and spin velocities are the same as in 
Eq.dini) and 



K. 



1+3411 ±94±- 92\\ T92± 
1 + 54|| ± 9i± + 92\\ ± 52_L - 2gi| 



1/2 



(79) 



C. Second order 

We next demonstrate how the renormalization of the 
occurs in the bosonic language, and how the "univer- 
sal" term in S with gfj_ emerges within the sine- Gordon 
model. Expanding Eg. ((70)) to order Sf_^^ and performing 
the averaging, we obtain the second-order piece in S 



E^'^^ = - J dxi j dx2 j (ix3Vi3l^23 cos[2fci^ {xi - X2] Jt{xi,X2)^ 



(80) 



where and all spatial integrals are cut at small distances by a. 

Again, this is nothing more that the two-bubble diagram 
2b) in FiglH written in the x, t space. The integration 

Jr{x,x')^\ dT']l2kF {x,T)Il2kp {x' ,t) , (81) 
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over T is readily performed 

coth [ttT ( 



4 sinh(27rr|a;|)sinh(27rr|a;'| 



There are two contributions to the term in S^^^ : one 
comes from large distances |a;| ~ \x'\ ~ and another 
one comes from distances of the order the interaction 
range. For the first contribution, the requirement that 
the potential must have a finite range is irrelevant, and 
the interaction in (|80p can be safely replaced by the delta- 
function V (x) = 2Trgi±S {x). We then obtain 



7(2) _ 



dx 



coth {2ttTx) 
sinh^ {2nTx)' 



J.2 



2^^-^sinh2 (27rTa)' 
Expanding the last result for Ta ^ 1, we obtain 



^(2) _ 



const 



(82) 



(83) 



This contribution is of the same magnitude but opposite 
in sign to the cutoff-independent part of 2 in Eq. (|T7)) . 
As the second contribution is expected to come from dis- 
tances smaller than T^^, we expand Eq. ((5T|) for T ^ 
and keep only the T dependent term 



Jr (x, x') 



J.2 



i\x\-\x'\f 



A8tt \x\ \x'\ {\x\ + \x'\)' 



Introducing new variables ^ = x — x' , rj — {x + x') /2 and 
performing elementary integrations, the resulting contri- 
bution to S can be represented as a sum of two terms 



-(2) 



— / dr^Wi?^) COS i2kFv)F+ in), 

rp2 roo 

— / d^W iO cos i2kF0 F- iO , (84) 

l^TT Jo 



where 



F+ iv) = log 



-2 + 4^, F_(0 = log 



(e/2 + a) 



(85) 

and W{x) = J dyV{x + y)V{y). The universal, cutoff- 

independent part of 2^ comes from the constant term 
(—2) in the function (rj) . It is of the opposite sign and 
twice larger than the contribution in Eq. (|83p . Combin- 
ing these two contributions together, we obtain for the 
universal part of 2 



^(2) _ 



The remainder of ^ is a cutoff-dependent part. To 
calculate this part, we consider two model interactions. 



The first one is consistent with the assumption used in 
the previous Sections (and also in g -ology, in general) 
that the backscattering part of the interaction is peaked 
near 2kF, i.e., the interaction oscillates in real space with 
period ir/kp- A model which describes this behavior is 

V (x) = 9iJ-^r~T^ {2kFx) . 
x^ -f 0^ 

The scale h equals to the bosonic cutoff A^^ introduced 
earlier. The assumption Ab <SC Af corresponds to the 
condition 6 3> a. Expanding functions F± for ^, ^ a 
and neglecting the exponentially small terms (of order 
exp {—2kpb)) as well as terms proportional to powers of 
a, we arrive at 



77(2) 



4r2 



nonuniv 



e 



^ log — 

f {2bf 2a 



ttT^ b 
^5i±log-, 
6 a 



where we used that dxlogx/ (l -I- x^) — 0. Combin- 
ing the universal and non- universal parts together, we 
obtain 



^(2) ^ _1 



l-21og- 
a 



(86) 



which coincides with Eq. (|T7)l upon identifying log - = 
ib. 

We see that the logarithmic renormalization of the 
backscattering coupling is reproduced within the sine- 
Gordon model. However, this result could have not 
be obtained for a local interaction. The interaction 
must have a finite range, which is larger than the short- 
distance cutoff in the theory. 

Another model, which we consider for completeness, 
corresponds to a long-range potential, i.e., to an inter- 
action peaked near q = in the momentum space. To 
describe this situation, we choose 

V{x)^^ 



TT x'^ + b'^ ' 



and assume that b ^ a 



so that the 2kF compo- 



nent of the potential V{2kF) = u exp(—2kpb) is exponen- 
tially small. Such interaction is not considered in the g 
-ology, and we will not express its parameters in terms of 
g -couplings. If backscattering is neglected completely, 
the problem is exactly soluble either via bosonization 
or Dzyaloshinskii-Larkin diagrammatic formalism^. It 
turns out that, somewhat surprisingly, corrections to the 
exact solution are small not exponentially but only alge- 
braically, in parameter 1/kpb. The reason is that loga- 
rithmic terms in functions F± [cf.Eq. ([55|) ]. which reflect 
correlations in motion of free fermions, introduce branch 
cuts into the integrals. The contribution of these branch 
cuts to the result is much larger than the exponentially 
small contribution of the poles in the interaction poten- 
tial. Evaluating the integrals and keeping only the lead- 
ing terms, we arrive at 



^(2) 



487r2' 



sin4fci?a ^ 2eb 
bkp a 



cos Akpo, 
2kj,ab 



O 



1 



kia%^ 
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where e = 2.718... The universal term, which is propor- 
tional to the 2k F component of the interaction is expo- 
nentially small and we do not retain it here. 

In the opposite case of a short-range interaction, i.e., 
for 6 <C a ~ /c 
term 



^(2) ig given entirely by the universal 



77(2) 



I2n 



IV. CONCLUSIONS 

In conclusion, we performed a detailed analysis of the 
temperature dependence of the specific heat for a ID in- 
teracting Fermi system. We used the g— ology model, and 
carried out a perturbative expansion in the couplings in 
the fermionic language. We have shown that, to first two 
orders in the interactions, the specific heat is expressed in 
terms of the non-running couplings in the RG -sense. The 
(74, and the 52 vertices appearing in C(T) are just bare 
vertices, while the backscattering gi vertex is the effective 
one, renormalized by fermions with momenta between 
fermionic and bosonic cutoffs. The running backscatter- 
ing amplitude on the scale of T appears in the specific 
heat only at third order in perturbation theory. The log T 
renormalization of the specific heat at the lowest T, ex- 
pected from the RG flow of the coupling constants, then 
only occurs at the fourth order in the perturbation the- 
ory, and the T-dependence of the specific heat follows the 
RG flow of the cube of the backscattering amplitude, in 



agreement with previous studies. We explicitly demon- 
strated that the absence of the logarithmic corrections 
below fourth order is due to cancellation of log T terms 
coming from low energies, of order T, and high energies, 
of order of the ultraviolet cutoffs in the theory. We also 
showed how the diagrammatic results can be obtained 
within the sine-Gordon model. 
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V. APPENDIX A 

In this Appendix, we derive Eq. for 
X = / dgllg;.^ (g, fi). The polarization operator 

n2fc^((7, ri) is given by (^5]) . It is convenient to split X 
into three terms X ^ Xi + X2 + X-^ as 



J 



Xo 



dq log 



4A? ' 



4A^/ .0 



dxdkknp{k) 



1 



1 



[q-inf-Ak^ [q + iO.Y-Ak'^ 
n 2 



(87) 



As we are only interested in a finite T contribution, we Using 
can safely subtract log^ from Xi. The 



rest is ultraviolet-convergent, and we can integrate ex- 
plicitly over q by setting the upper limit of the integral 
to infinity. We obtain 



T^|r!|log|f]| 



■logT 



TT 



77- ' * 



log 2 



we obtain 



X log 2x 
sinh'^ X 



X, = 



2T' 



0.5803, 



T 



(89) 



(90) 
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where i? = | — log 2 — -^h = 0.454. The second term, 
X2, contains contributions both from small q, of order T, 
and from large q, of order Ab. It is convenient to split the 
momentum integral /p ^ into — . The first integral 
can be easily converted into the integral over the whole 
real q axis. The poles in q at any finite O are located in 



the same half-plane, and the g— integral is nonzero only 
because of the branch cut in the logarithm. Choosing 
the integration contour as shown in Fig.[7l evaluating the 
momentum integral, performing the frequency sum, and 
adding a separate contribution from J7 = 0, we obtain 



r2 



f°° £7^ + q"^ f°° f 1 1 \ 2r^ 



The integral over large g > Ab involves also large fre- 
quencies ft q, and the frequency sum can be safely 
replaced by the integral. Typical fermionic momenta, fc, 
are of order T and, therefore, much smaller than q and 



log^-S-(0.5 + log2) 

. m 



fl. Neglecting k in the denominators of the integrand, 
and performing three independent integrations (over fc, 
il, and q), we obtain 



dq log 



dkknp{k) 



1 



1 



(q-iny^ 



8 f 1 2T'^ 

- dkknAk) J dn j^^ ^'^(^3-^^)1 - - — 



-4fc2 

, 2Af 
'°^A^ 



{q + iny 



4fc2 



-(2+ log 2) 



(92) 



Combining the two contributions, we obtain 



2r2 



log 



T 

2AT 



log 



2Af 
AT 



B 



(93) 



In the third term, X^, the 2D integral over q and fl is 
ultraviolet convergent, and we can safely set the upper 
limit of the q— integral to infinity. The integral again 
has separate contributions from 7^ 0, and from = 0. 



The contribution to X^ from finite frequencies is eval- 
uated straightforwardly by closing the contour of the 
q— integral in the upper or lower half-plane. We obtain 
(r2/3)(7 - 10 log 2). The evaluation of the contribution 
from = requires special care because of the poles 
which are avoided by replacing f2 by i5. The correspond- 
ing contribution to S3 becomes 



—^tJ dq j dkknp{k) j dppnp{p) 
= 4T ^ n| {x)dx = 4T2 ^log 2 - i 



1 



1 



{q-iSy-Ak"^ [q + iSy - Ap^ J \{q ~ iSy - Ap^ [q + iSy - Ap' 



1 



1 



r 



(94) 



We emphasize that the integral in the r.h.s. of (|M|) comes Collecting ([TO]) . and (I^S)) . we obtain 

from an infinitesimally small region where |fc — p| ~ 5. 



Combining the two contributions to X^, we obtain 



1\ 

X3=(log2+- — . 



(95) 



Eq. dM]) coincides with Eq. 



(96) 
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-i^l-2k • •-i^2+2k 



FIG. 7: Integration contour for Eg. (191) . 



VI. APPENDIX B 



In this Appendix, we derive the resuh for Y = 
^ J2n I dl^2kF ^) logarithmic accuracy. We as- 
sume that T is small, such that T <C Ab,Af and that 



Ab ^ At, and collect terms logarithmic in T/Af, and in 
Ab/Af. The computational steps are the same as in Ap- 
pendix A: we use the fact that Il2kF given in is the 
sum of two terms and split Y into Yi + Y2 + Y^ + Y4, 
where 



y,=ty: 



Y. 



■log 



2 I „2 /-oo 



2 + 



4A^ .0 



f 

dknp{k)k 



\{q~iny~4:k^ {q + iVlf ~ Ak"^ J ' 



Yi ~ 3"-^^^ / d,q dxnpi 



{x)x 



1 



1 



(97) 
(98) 
(99) 
(100) 



One can easily make sure last term I4 is non-logarithmic 
and can be neglected. 

The momentum integral in Yi is infrared divergent. 
However, we only need the thermal part of Yi. To ex- 
tract it, we subtract from the integrand in Yi its value at 
$1 = 0, i.e., log -fi^- This makes the momentum integral 
finite. Evaluating it and then performing the summation 
over frequency, we obtain 



log^ ^- 0.909 log £ 



(101) 



where dots stand for 0{T^) terms. The number, 0.909, 
as well as other numbers below are expressed in terms of 
convergent ID integrals. 

In Y2, the cutoff in the integration over qi is relevant. 
Splitting the g— integral into /q into — and eval- 
uating each of the two terms separately in the same way 



as in Appendix A, we obtain 



2 r 



log' =^ - 3-2!l5 log + log' ^ - log 

ZAf Af Ab Ab 



(102) 

The result from Y3 can be readily obtained from the ex- 
pression for X3 in Appendix A, as to logarithmic accu- 
racy we can replace log ^Ip^ in ([99]) by 2 log (T/Af). We 
then obtain 



y2 y 

Y3 = 2.386 — log—. 

TT Af 



(103) 



Combining pOTj) . ^02]) . and p03)) . we obtain that all 
log (T/Af) terms are cancelled out, and 



TT 



, 2 Af Af ' 



(104) 



20 



Eq. (fT04)) coincides with ((53)l . 

Another backscattering diagram which could possibly 
give rise to logarithmic terms is diagram 3b) in Fig. [5l 
For a local interaction, it reduces to a cube of the Cooper 
bubble, which in ID coincides with n2fcj^ up to the overall 
sign. However, one can easily verify that for Ab ^ Af, 
this diagrams does not contain log(Af/Ab) terms. In- 



deed, the cutoff induced by the interaction imposes the 
restriction on three out of four momenta and frequencies 
in the fermionic lines. The 2D integral over the remain- 
ing momentum and frequency involves all six fermionic 
propagators and is confined to the lower limit. This im- 
plies that all variables are of the same order, and there 
is no space for a logarithm. 
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